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1.   INTRODUCTION 

As  the  requirements  for  computer  circuit  performance 
becomes  mere  demanding,  there  is  an  ever  increasing  need  for 
a  more  efficient  method  of  describing  transistor  circuits. 
Before  this  need  can  be  met,  a  transistor  model  must  be 
available  which  will  permit  accurate  prediction. 

There  are  numerous  models  available  which  (3)  permit 
accurate  prediction.   Many  parameters  associated  with  these 
models  are  not  easily  evaluated,  and  the  evaluations  often 
require  facilities  and  measurement  techniques  which  are  not 
commonly  available.   These  models  cannot  be  easily  simulated 
on  the  computer;  such  a  simulation  often  requires  advanced 
rype  computer,  and  program  execution  times  are  lengthy.   In 
short,  these  models  are  cumbersome  to  use  in  practice.   There 
is  need  for  a  model  which  is  based  on  theoretical  principles 
but  which  can  be  easily  used  in  practice.   Such  a  model  1) 
should  be  a  good  compromise  between  accuracy  and  tractibility, 

2)  should  have  the  property  of  ease  of  evaluation  of  its  para- 
meters with  conventional  facilities  and  reasonable  precision, 

3)  should  permit  simulation  on  the  computer  easily  with  a 
reasonable  execution  time,  4)  should  be  valid  for  large  as 


well  as  small  imput  signals.   Large  signal  input  can  be  as 
large  as  to  derive  the  transistor  into  saturation,  but  should 
not  be  so  large  to  cause  high  injection  effects. 

The  model  to  be  described  in  this  report  is  not  a  complex 
model.   In  fact,  it  has  been  around  for  quite  some  time. 
However,  with  the  addition  of  the  nonlinear  characteristics 
and  with  careful  device  characterization,  this  model  can 
provide  the  necessary  means  of  predicting  response  to  large 
as  well  as  small  input  signal. 

The  characteristics  of  the  p-n  junction  have  been  studied 
by  many  authors  (2).   By  concentrating  study  on  only  one  junc- 
tion, one  can  understand  more  readily  some  of  the  more  impor- 
tant transistor  characteristics.   In  developing  this  model, 
each  junction  of  a  common  base  transistor  is  studied,  while 
the  other  terminal  is  shorted  to  base.   The  use  of  superposition 
is  made  in  obtaining  the  model.   Because  the  model  is  nonlinear, 
the  use  of  superposition  may  at  first  be  considered  invalid. 
The  model  is  to  be  simulated  on  the  computer,  and  the  computer 
obtains  the  solution  to  the  differential  equations  by  means 
of  numerical  integration.   By  choosing  small  enough  intervals 
of  integration,  nonlinearity  can  be  observed  as  linearity 
within  each  interval.'  The  program  re-evaluates  the  nonlinear 
elements  for  each  given  interval.   Then  using  the  system  of 


nonlinear  differential  equations  describing  the  model,  in 
conjunction  with  the  internal  subprograms,  the  value  of  the 
voltages  corresponding  to  that  interval  are  predicted  and 
corrected.   If  the  predicted  value  is  not  within  specified 
allowable  error  from  the  corrected  value,  then  the  program 
modifies  the  interval  of  integration  and  repeats  the  cycle 
until  the  specified  allowable  error  is  met.   Thus  use  of 
superposition  is  validated. 

Most  of  the  methods  described  in  this  report  to  evaluate 
the  parameters  of  the  model  are  standard  methods;  however, 
they  are  modified  when  necessary. 

Although  the  model  is  tested  with  fast  rise-time  input 
pulses,  the  model  has  more  general  applications.   The  input 
pulse  tests  are  more  severe  than  simple  a-c  or  d-c  input 
tests.   By  fourier-series  expansion  it  can  be  shown  that  the 
pulse  has  both  a-c  and  d-c  components. 


2.   EQUIVALENT  CIRCUIT 

The  model  being  derived  is  based  on  the  Ebers  and  Moll's 
equations  for  transistors  (1) ,  Equation  (1)  and  (2) .   The 
direction  of  the  currents  and  polarity  of  the  voltage  terminals 
are  shown  in  Fig.  (l) . 
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Fig.  1.   The  direction  of  current  and  polarity  of  the 

voltage  which  corresponds  to  Equations  (l)  and  (2) 


Where : 


I_  =  Emitter  Current 

E 
Ir   -   Collector  Current 

VE  =  Emitter-to-Base  Voltage 

Vr   =  Collector-to-Base  Voltage 


IE   -   a 


H 


(eqVEBAT   _    1}    .    a       ,eqVCBAT    _    1} 
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I„   =   a 
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(eqVEBAT    _    1}     _    a       (eqVCBAT    _    l} 
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(1) 
(2) 


Where : 


q  Electron  Charge 

k  Boltzmann's  Constant 

T  Temperature 

a,-.,  a29/  a2i*  and  a12  are  Constants 


The  constants  a   ,  a,,,,  a--,  and  a22  are  functions  of 
transistor  material  and  junction  temperature.   The  meanings 
given  to  these  constants  in  this  paper  are  different  from  the 
meanings  that  are  given  in  the  reference  (1) .   This  will  be 
further  discussed.   Equations  (1)  and  (2)  are  derived  solely 
from  the  diffusion  equation  for  the  d-c  condition.   The 
assumptions  made  in  deriving  Equations  (1)  and  (2)  are:   that 
drift  current  is  negligible  compared  to  diffusion  current, 
and  that  there  are  no  mobile  charges  in  the  depletion  layers 
(7) .   The  above  assumptions  are  valid  for  the  first  derived 
model.   The  model  obtained  in  this  manner  will  be  called  the 
intrinsic  model.   The  intrinsic  model  will  be  further  com- 
pleted by  introducing  some  of  the  extrinsic  elements.   The 
extrensic  elements  considered  are:   depletion  junction 
capacitance,  bulk  resistance,  and  header  capacitances.   The 
model  resulted  by  introducing  extrinsic  elements  into 
intrinsic  model  is  called  extrinsic  model.   The  extrinsic 
model  further  be  completed  by  introducing  additional  non- 
linear elements  to  the  model  to  account  for  base  width  modu- 
lation and  other  effects  caused  by  forward  baising  collector 
junction,  and  this  model  will  be  called  complete  model. 

The  solution  to  Equation  (1)  and  (2)  for  I„  and  Ic  are 
obtained  by  superposition.   This  is  done  1)  by  setting  the 


voltage  VCB=  0  and  obtaining  the  IE  and  Ic  only  due  to  VgB; 
2)  by  setting  the  voltage  VEB-  u  and  obtaining  the  IE  and  Ic 
only  due  to  VCE;  3)  the  IE  and  Ic  that  satisfies  equations 
(l)  and  (2)  are  the  summation  of  IE ' s  obtained  from  condition 
(1)  and  (2)  and  summation  of  Ic ' s  obtained  from  condition  (1) 
and  (2) . 

The  value  of  V"cs  is  set  equal  to  zero  by  shorting  the 
collector  to  base.   The  equations  (1)  and  (2)  are  then  reduced 
to  equations  (3)  and  (4)  respectively. 


IE1   =   au(e^BAT   _    1}  (3) 

IC1   =   a21(e1%AT   _    d  (4) 


Equation    (3)    is    analogous    to    the    voltage    v.s.    current   charac- 
teristics   of   a  p-n   junction   diode.      The    current-voltage 
characteristics    of   a  p-n  junction   can  be    represented  by   equa- 
tion   (5)  . 


I   =    Is    (e^T   -    1)  (5) 

where : 

Is  back  saturation  current  of  a  p-n  junction, 


Is  is  function  of  cross  sectional  junction  area,  diffusion 
constant,  diffusion  length  and  dopping  of  acceptors  and  holes 
in  p  region  and  n  region. 

A  close  comparison  between  equations  (3)  and  (5)  yields 
that  constant  a,  j_  is  some  sort  of  back  saturation  current  of 
a  p-n  junction.   So  a]_j_  will  be  defined  as  Ise*  tne  back 
saturation  current  of  the  eroitter-to-base  junction  when  the 
collector  is  shorted  to  base.   By  substituting  ISE  in  equa- 
tion (3) ,  equation  (6)  is  obtained. 


*E1    =    *SE     (e^VEBAT    _    l}  (6) 

By  Kirchboff's    current    law: 

IE    =   -(IC   *    1b)  (7) 


By  dividing  equation  (4)  into  equation  (6)  a2]_  is  obtained, 
equation  (8) . 


a21  =  ^L.  ISE  (8) 

-"■El 

A±    is    defined  as  ^£1     =  A-,  (9) 

JE1 


Further  manipulation  results  equation  (10) 


IC1   =  A1ISE    (eqVEBAT   "    1)  (10) 


0r  JC1    =    A^El 


Equations  (6)  and  (11)  suggest  that  the  below  model  can  be 
derived,  Fig.  (2) .   The  model  Fig.  (2)  suggests  that  for 
understanding  the  behavior  or  the  emitter  junction,  one  can 
investigate  the  modified  diode  of  Fig.  (2) .   Hence,  the 
characteristics  of  the  diode  will  be  investigated.   The 
voltage  vs.  current  characteristics  of  the  diode  p-n  junction 
is  shown  in  Fig.  (3),  and  for  an  ideal  diode,  the  functional 
relationship  representing  this  V-I  characteristic  is  given 
by  equation  (12) . 

The  functional  relationship  describing  the  conductance 
of  the  junction  due  to  the  diffusion  mechanism  is  defined  in 
equation  (13) . 


gd  =  I  =  ^(y  i)  d3) 
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V  is  applied  Voltage. 
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Fig.  2.   Model  representation  of  emitter  junction  when 
VCB  "  °- 


Fig.  3.   V-I  Characteristic  of  Diode 


I  =  Is(e  V  -  1) 
=  q/kT 


(12) 


10 


"El 


V 


EB 


B 


•ci 


k»  ^Cdb  Qlx 


VCB    =    ° 


■El 


-5   B 


Fig.    4.      Model   representation   of   emitter   junction  with 
nonlinear   elements,    when  VCB   =    0. 


Where : 


G         ■    W*»VEB    -    1) 
DB  VT 


(15) 


EB 


On«   =  TV    ISKe*VEB 


'DB    ~    '  b    ■LSE 


(16) 
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The  functional  relationship  describing  the  diffusion  capaci- 
tance of  the  junction  is  given  by  equation  (14). 

CD  "  lbW   =  'b  Ise&V  (14) 

"T^   the  effective  minority  lifetime  in  the  base. 


The  model  in  Fig.  (2)  can  be  presented  as  shown  in  Fig.  (4) 

Now,  by  setting  the  value  of  Vj?b  equal  to  zero  in 
equations  (1)  and  (2) ,  characteristics  of  the  collector 
junction  may  be  obtained  which  are  similar  to  the  emitter 
characteristics  when  VCg-  0. 


*E2  =  -*12(^VCB  ~  1) 


(17) 


IC2  "  -a22(e^VCB  _  1}  (1Q) 

t*   =  qAT 

A2    is    defined  as  A2   =    IS2  (19) 

JC2 


A  close  comparison  hetween  equations  (5)  and  (la) 
yields  that  constant  a22  1S  some  sort  of  back  saturation 
current  of  a  p-n  junction.   Then  a22  will  be  defined  as  Igc 
the  back  saturation  current  of  the  collector  to  base  junction 
when  the  emitter  is  shorted  to  the  base.   By  substituting 


12 


Isc  in  equation  (18) ,  equation  (20)  is  obtained. 


*C2  =  -Isc(^VcB  -  1)  (20) 


Further  manipulations  similar  to  the  manipulation  done  for 
the  emitter  junction  yield  equations  (21-22)  for  collector 
j  unction. 


JC2  =  -ZSC^03   ~    1)  (2D 

*E2  =  A2  *C2  <22> 


The  equations  (21)  and  (22)  suggest  that  model  Fig.  (5)  can 
be  derived.   The  model  in  Fig.  (5)  can  be  presented  as  shown 
in  Fig.  (6)  by  introducing  nonlinear  elements. 

Now  the  solution  to  the  equations  (1)  and  (2)  are: 

*E  =  -^El  +  IE2  (23) 

Ic  =  ICi  4  IC2  (24) 


J22 
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Fig.    5.      Model   representation   of  collector   junction 
when  VEB   =    0. 
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Fig.  6.   Model  representation  of  collector  junction 
with  nonlinear  elements,  when  VEB  =  0. 
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Therefore : 

*E  =  ^SE(^VeB  "  1)    '  (25) 

Ic  =  ^1IE1   -  Isc(e^VcB  -  1)  (26) 


Equations  (25)  and  (26)  suggest  that  model  in  Fig.  (7)  can 
be  derived.   The  Model  in  Fig.  (7)  can  be  represented  as  in 
Fig.  (8)  by  introducing  nonlinear  elements. 

An  equivalent  circuit  for  the  common  emitter  configura- 
tion of  the  equivalent  circuit  of  Fig.  (8)  will  be  obtained. 
The  first  step  is  to  redraw  the  circuit  of  Fig.  (8)  into  a 
"T"  equivalent  circuit  as  shown  in  Fig.  (9) .   Further  manipu- 
lations of  Fig.  (9)  suggest  the  below  equations. 


IB  -  IE1  -J-  A2IC2  "  A1*E1  ~  JC2  *    (35> 

Or    IB  =  -IE1(1  +  A-l)  -  Ic2(l  +  A2)  (36) 

And   Ic  =  A1IE1  -f  Ic2  (37) 

Or    Ic  =  Ic2(l  ~  A2)  +  ^C2A2  +  A1IE1  ^38^ 


Equations  (36)  and  (37)  suggest  Fig.  (10) . 

A-^  and  A2  can  be  related  to  hFE  and  hFEI  respectively. 
This  is  done  as  follows.   The  transistor  is  operated  inversely 
and  VEB  is  set  equal  to  zero. 
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Fig.  7.   Model  representation  of  a  common  base 
transistor. 


IE1  =  ISE(eWEB  -  !) 


IC2  -Isc(e^CB  -  1) 


IE2  =  A2*C2 
IC1  =  ^El 


(27) 
(28) 
(29) 
(30) 
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Fig.  8, 


Intrinsic  Model  representation  of  a  Common  base 
configuration  of  transistor  with  nonlinear 
elements. 


Where : 


GDB  '  -j2ik" 


V 


EB 


CDB  =  'b  XSEe 


XVEB 


X =   g/kT 

IQr(e^VCB  -  1) 


-  XSC 


SB 


V 


CB 


CSB  =  rbhFE  ISCe/<VCB 


f   =  qAT 


hFE  -  short-circuit  current  gain  for  common-emitter 
mode  of  transistor  operation. 
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Fig.    9.       "T"    Equivalent  Circuit   of   Fig.    8. 
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Fig.    10.       Intrinsic   Common  Emitter   Equivalent   Circuit 
of  Transistor. 


*E1    =    W^ BE    "    « 


IC2    =    -Iscte^CE    "    VBE)-    1) 


(39) 
(40) 


The  parameter  hFEI  is  defined  as  short-circuit  common-emitter 
current  gain  for  a  transistor  operated  inversely. 


'M   =  hFEI 
•B 


(41) 
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A?    =  £e  (42) 

Or        A      =    -   ff  (43) 
E          XB 

0r         A2    =    -  TTH5  -  lThFEI                                           l44) 


Now  VCE  is  set  equal  to  zero,  therefore 


Ir  Al  /   v 

-°-  =  — ± =  hFE                         (45) 

IB  1  ■-  Ax 

Or    An  =  ^FE  ,  (46) 

1  I  +  hFE  V   ' 


The  Intrinsic  Model  representation  of  a  common-emitter  con- 
figuration of  transistor  with  nonlinear  elements  is  shown  in 
Fig.  (11) . 

The  model  in  Fig.  (11)  is  based  on  the  diffusion  mechanism 
only,  and  is  called  an  intrinsic  model.   This  model  can  be 
further  completed  by  introducing  extrinsic  elements  in  the 
model  which  are:   junction  depletion  capacitance,  bulk  resis- 
tances, and  header  capacitance.   Although  the  extrinsic  ele- 
ments are  not  derived  from  diffusion  equation,  they  may  be 
dependent  on  the  junction  current  or  junction  voltage.   This 
will  be  discussed. 
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Fig.    11.       Intrinsic  Model    for   a   Common-Emitter 
Configuration. 


Where : 


G^   = 


_    ISE(e-*VBE   "    1) 


D  (hFE    4    1)VBE 


-?v  I< 


-*V 


CD  YhFS    1    1) 


BE 


iSc^nVBE~  Vce)  - 1) 

GS      ■     (hFEI    +    1) (VCE    -    VBE) 


CS    = 


_rbhFE    e~nV^    ~   V^] 


(hFEI    +    1) 


i!    =   hFE   VBEGD 


I2    =    hFEI   GS(VBE    -   VCE) 
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The  variation  of  the  depletion  layer  capacitance  of  the 
diode  p-n  junction  as  a  function  of  the  applied  voltage 
depends  on  the  nature  of  the  charge  distributed  on  the  junc- 
tion.  This  capacitance  for  the  abrupt  charge  density  distri- 
bution is  in  the  form  of  equation  (47) . 


CJ  =  w  -  v>*  C47) 

Where :   C^  is  a  constant  function  of  the  dopping  level, 
area  of  the  junction  and  permittivity  of  the  junction. 
Vp    is  the  built-in  depletion  layer  voltage  which  is  a 
function  of  serai-conductor  material.   V  is  the  externally 
applied  voltage  across  the  junction. 


The  depletion  layer  capacitances  for  the  linearly  graded 
charge  density  is  in  the  form  of  equation  (48) . 


Cl 

CV^  -  Vjl/3 


Cj  =  ,„   ^t.  „  (48) 


The  charge  distributions  used  in  the  transistor  junctions  are 
usually  either  abrupt  or  linearly  graded.   However,  they  are 
not  ideally  abrupt  or  linearly  graded,  so  that  the  value  m  is 
chosen  instead  of  the  powers  1/2  or  1/3  for  (V^-  V). 
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Thus  CT  =  El _ 

J    (V^  -  V)m 


It  should  be  pointed  out  that  the  assumption  was  made  that 
there  are  no  mobile  carriers  in  the  depletion  layer. 

The  bulk  resistance  of  the  emitter  region  or  the  collec- 
tor region  is  obtained  from  the  conductivity  of  the  emitter 
region  or  the  collector  region  respectively.   However,  these 
resistances  are  small  enough  to  be  ignored  for  all  practical 
purposes.   The  bulk  resistances  of  the  base  region  is  not  only 
due  to  the  conductivity  of  the  base  region,  but  also  depends 
on  the  recombination  of  the  carriers  in  the  base  region 
which  is  strongly  dependent  on  the  magnitude  of  emitter 
current.   This  phenomena  is  not  well  understood. 

The  header  capacitances  are  due  to  the  lead  wires  to 
transistor  and  the  transistor  encopsilation.   These  capaci- 
tances could  affect  the  performance  of  the  transistor  at  high 
speed  application. 

The  model  including  the  extrinsic  parameter  for  the 
common  emitter  configuration  is  shown  in  Fig.  (12) . 

The  functional  relationship  representing  the  parameters 
of  Fig.  (8)  for  PNP  and  NPN  transistor  are  on  page     of  this 
report. 
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Fig.    12.      Extrinsic  Model    for  Common  Emitter  Configuration. 
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The  parameters  corresponding  to  Fig.  (9)  for  PNP  and 
NPN  transistors  are  given  below: 

PNP  NPN 


—       ^El  —       ^El 

°E  ~  (V  E  -  VB,E)ra  CE  "  (VB.E  -  V  E)m 

c  =         ^ci  c  =         ^ci 

°    (V  c  -  (VCE  -  VB,E))n  C    ((VCE  -  VB.E)  -  V  c)n 

oj  q-XVr'f  t         *vB'E 

c      _   /fa    ISE    e    qvB   E  c  _      b    ISEe 

D     (hFE  4  1)  D     (hFE  4.  i) 

t   /  ~VB'E    i\  T   /  &VB'E    i\ 

r      =  ISE(e      ~  1)  r,  =  ISE(e      -  1) 


(hFE  4  1)VB.E  (hFE  4  1)VB.E 

.rischFE  (e-^CE-VB'E)  -  1)       .TbIsrhFR  (e^^'^  -  1) 
S  -  b         (hFEI  4-D  CS  (hFEI  4  D 

Qs  =  Isc(e^tVCE-VB-E)  -  l)         Gg  =  Isc(e^VCE-VB'E)  .  1} 


(hFEI  +  1)  (VB.E  -  VCE)  (hFEI  4-1)  (V    -  VB,E) 


CE        Emitter  junction  depletion  capacitance. 

CE1       The  value  of  CE  when   (VE  -  VBiE)   =1. 

m         The  parameter  related  to  the  CE,  evaluated  emperically. 

Cq,    Cq-^i    n  Collector  parameters  which  are  analogous  to  the  CE, 

CE^,  and  m  respectively. 
CD        Emitter  junction  diffusion  capacitance. 
T-b  The  effective  minority  carrier  lifetime  in  the  base. 

It  is  a  function  of  q/kT  for  emitter  junction. 
ISE   Back  saturation  current  for  emitter  junction. 
hFE   Short-circuit  current  gain  for  common-emitter  mode,  of 

normal  transistor  operation.   Equation  (45) . 
vyE   The  barrier  potential  of  the  emitter  junction. 
t1 '    ■"■sc  V^C   Collector  parameters,  which  are  analogous  to  the 

ft  ,  IsE'  and  vj#E'  respectively. 
hFEj      Short-circuit  common-emitter  current  gain  for  a 

transistor  operated  inversely.   Equation  (41) . 
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3.   DEVICE  CHARACTERIZATION 

The  techniques  used  to  characterize  the  device  will  be 
discussed  here.   It  will  become  apparent  later  from  the  sensi- 
tivity studies  that  not  all  device  parameters  discussed  will 
need  to  be  measured  on  each  and  every  device  (5) .   In  practice 
nominal  values  will  suffice  for  a  given  family  of  devices. 

Because  of  the  similarity  in  which  the  emitter  and  collec- 
tor junctions  are  represented  in  the  model  it  is  sufficient  to 
describe  only  the  measurements  performed  on  the  emitter  junc- 
tion in  detail.   The  few  differences  that  exist  in  the  measure- 
ment techniques  for  the  collector  are  appropriately  cited. 

The  parameter  I__  and  &  are  determined  with  the  emitter 

biased  in  the  forward  direction  and  the  collector  shorted  to 

base.   The  emitter  input  is  a  ramp  input.   Although  the  value 

of  q/kT  is  attached  to  the  X  iR  most  device  treatises,  it  is 

found  experimentally  that  the  X  is  equal  to  the  .  SL  where  a 

1  u  akT       - 

is  experimentally  found.   The  existence  of  a.  can  be  explained 
by  the  approximation  used  in  evaluating  GD  and  C  .   Further 
manipulation  of  equation  (5)  reveals  equation  (50) . 
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*E  =  ISe(^V'-  1)  (5) 

°r    IB  =  ISE^V       e*V»l  (50) 

And   logxolE  =  lo<3lQ1SE   +  v(^lO9i0e^  ^50^ 


The  value  of  ISE  then,  is  obtained  by  plotting  log,QIE  vs.  V. 
The  intersection  of  the  extrapolated  line  with  the  log-mlg 
axis  is  the  value  of  log  ISE  and  the  slope  of  the  line  is 
^log-jQe.   A  plot  of  the  V"EB  vs.  IE  is  shown  in  Fig.  14.   It 
is  incorrect  to  use  the  terminal  voltage  VEB  for  the  value  of 
V  in  equation  (50) ,  because  the  V  in  equation  (50)  is  the 
intrinsic  voltage  applied  across  the  junction.   The  value  V 
in  terms  of  V  „  is  given  in  equation  (51) . 


Where  RSE  is  emitter  bulk  resistance. 


The  modified  V-I  characteristic  is  shown  in  Fig.  (13) .   The 
graphical  technique  used  to  evaluate  y  and  ISE  for  the  transis- 
tor X  is  shown  in  Fig.  (14) .   (See  Appendix  1  for  information 
on  Transistor  X) . 

The  Barrier  junction  potential  M$   can  be  evaluated  by 
equation  (52) . 
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V  (volts) 


Fig.  13.   The  Diode  V-I  Characteristic. 


Solid  line  is  terminal  voltage  v.s.  current. 

Dash  line  is  corrected  value  of  voltage  v.s.  current. 
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Fig.  14.   Graphical  technique  to  evaluate  IqF  and  & . 
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Vjrf  =  H  in  £*£§  (52) 

.  q       Ni^ 

Na  =  density  of  acceptor  atoms. 

Nd  =  density  of  donor  atoms. 

Ni  =  carrier  density  in  an  intrinsic  semi-conductor 


Equation  (52)  is  very  cumbersome  to  use  in  practice.  Vtf   can 
be  evaluated  easily  by  approximation.   The  sensitivity  study 
has  revealed  that  the  V0   is  not  a  very  sensitive  parameter 
and  a  close  approximation  is  accurate  enough.   A  diode  at 
high  current  can  be  approximated  by  a  series  circuit  consis- 
ting of  a  battery  equivalent  to  barrier  potential,  a  resis- 
tance which  is  equivalent  to  the  diffusion  resistance  of  the 
diode,  and  an  ideal  diode.   The  extrapolation  of  the  slope 
to  the  V-I  characteristic  will  meet  the  V-axis  approximately 
at  Value  V0.      The  graphical  evaluation  of  this  barrier  poten- 
tial is  shown  in  Fig.  13.   The  nominal  value  of  .8  volts  can 
be  used  for  a  silicon  transistor  (5) . 

CE1  anc^  ra  are  obtained  from  capacity  measurements  on 
the  emitter  junction  (collector  open  circuited).   The  capaci- 
tance remaining,  after  the  contribution  of  the  lands  and 
header  are  removed,  versus  the  true  junction  voltage  is  plotted 
on  log-log  paper.   The  junction  voltage  is  the  applied  voltage 
less  the  contribution  of  the  barrier  potential.   Fig.  15  is  an 
example  of  how  the  values  of  C^i  and  m  are  obtained. 
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x  Experimentally  observed  point 


slope-  ra 


(VBE  -  0.8)  volts 


Fig.  15.   Graphical  technique  for  obtaining  the  value 


m  and  C  .. 
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A  computer  program  can  be  written  to  evaluate  C„n ,  ra  and 
Vo7  simultaneously.   This  can  be  done  by  sequential  least 
square  fit,  which  will  find  those  values  of  CE-,,  m  and  Vo7  to 
fit  experimental  data  and  equation  (49) . 

Short  circuit  common  emitter  current  gain  is  measured  in 
the  conventional  manner.  Fig.  16  is  an  example  of  the  depen- 
dency of  this  parameter  on  emitter  current. 

The  value  of  i-^   may  be  obtained  by 


T   =     1'22  i^\ 

'b    (1-%)  WN  ^53) 

Where  d,-^   is  low  frequency  of  short  circuit  current 

gain  for  common  base  circuit. 

Where  WN   (X  cut  off  frequency  for  the  common  base 

short-circuit  current  gain  for  a  normally  operated 

transistor. 


In  practice  it  is  cumbersome  to  use  equation  (53)  for  evalua- 
tion of  T   .      The  charge  method  technique  is  more  widely  used 
for  evaluation  /-,,    which  is: 


7^=  AQfi  (54) 


b   AIB 
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Fig.  16.   Common  emitter,  short  circuit  current  gain 
versus  emitter  current. 
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ff     is    defined  as    /      =  4_?B  (55) 


A^c 


By  manipulating  equations  (55)  and  (54) ,  equation  (56)  will 
result. 


£=■?!•  hFE  (56) 


B 


Commonly  /  is  found  first  and  then  V-q  is  calculated.  This 
is  done  because  the  error  involved  in  evaluating  /Q  is  less 
than  the  error  in  evaluating  /g. 

Base  spreading  resistance,  Rg  can  be  obtained  from  a 
small-signal  bridge  measurement  of  the  Rb'cc  Pro<3uct,  Fig.  17. 

The  values  of  Rfi  indicate  that  RR  is  a  function  of  the 
emitter  current  and  decreases  as  the  emitter  current  increases; 
however,  these  values  level  off  at  higher  emitter  current. 
For  the  transistor  X,  a  nominal  value  of  15  ohms  is  obtained 
for  the  emitter  current  beyond  3ma.   For  emitter  currents 
less  than  3ma,  Rs  increases  to  about  twice  the  nominal  value 
at  lma. 

ISC'///  CC1'  n  and  hFEI  are  determined  in  the  same  manner 
as  ISE/  X ,  CE1,  and  m  and  hFE  respectively,  with  the  exception 
that  the  emitter  and  collector  leads  are  interchanged. 
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enerator 


VE   Bias 


Fig.  17.   The  diagram  used  to  evaluate  Rfi, 


R  is  Known  Resistance. 
C^r  is  Known  Capacitance 
D   is  Detector. 


And: 


R, 


=  Rkck 
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4.   SIMULATION  ON  THE  COMPUTER 

The  model  in  Fig.  12  is  programmed  on  IBM  7090.   This 
program  can  be  used  on  the  IBM  709X  series.   The  program  is 
written  in  the  FORTRAN  II  language.   A  brief  discussion  of 
the  program  is  given  below.   (See  Appendix  B) 

The  time  response  of  the  circuit  is  found  by  forming 
node  equations  and  solving  for  derivatives  where  there  are 
reactive  elements  at  the  node.   This  yields  a  system  of  non- 
linear differential  equations.   These  equations  are: 


GBVB  "  GBVBP  +  <CBES  +  CCBS>VB  "  CCBSV0  =  0        (59) 

-VBGB.4  (GB  4  GD  4.  Gs)VBp  -.GSV0  *     (CE  4  CD  4  Cc  4  Cs) 
VEp  -  (Cc  4  Cs  4  CCBS)V0  =  0  (60) 

-GsVBp  4  GSV   -  CCBSVB  "  <CC  +  CS>VBP  +  (CCBS  4  Cc) 

VQ  =  -hFEVBpGD  4  hFEI  Gs (VQ  -  VBp)    (61) 


Where : 

VB  =  VBE 

VBP  =  ^E 

V0  =  VCE 

vb  "  m? 

^BP  =  dtVBP 

*o  =  ltvo 

Equations  59  -  61  are  solved  by  a  predictor-corrector  numerical 
integration  technique.   Nonlinear  elements  are  handled  by  com- 
puting their  value  at  each  interval  of  integration  step.   The 
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accuracy  and  running  time  of  this  program  is  highly  dependent 
on  the  initial  value  of  voltages,  the  integration  interval, 
and  the  maximum  allowable  relative  error  for  each  integration 
interval.   The  value  of  1.0E  -  7  has  yielded  accurate  results 
and  is  recommended  for  use.   The  interval  at  which  the  program 
prints  also  determines  the  minimum  integration  interval  which 
in  turn  affects  the  accuracy  of  the  solution.   Accurate  cal- 
culations of  the  initial  condition  of  the  voltages  (V  ,  VBp, 
and  VQ)  will  decrease  running  time  and  have  a  solutary  effect 
on  the  computed  response.   The  input  signal  to  the  transistor 
should  not  be  programmed  at  the  initial  time  Tq.      This  is  to 
allow  the  transients  introduced  by  inconsistent  initial  condi- 
tions to  dissipate  before  applying  the  input  signal. 

The  functional  relationships  describing  the  RB  and  hFE 
as  a  function  of  the  emitter  current  are  of  complex  nature; 
also,  these  relationships  are  not  very  representative  of  the 
phenomenon)  occurring.   The  values  of  RD  and  hFE  as  a  function 
of  the  emitter  current  are  given  by  the  table  of  values  to 
the  computer.   An  error  is  introduced  in  this  program  to 
reduce  the  computing  time.   This  is  done  by  evaluating  the 
present  value  of  RB  and  hFE  through  using  the  last  value  of 
IE.   This  error  is  of  no  consequence,  and  the  approximation 
is  justifiable,  because  the  variation  of  the  Rfi  and  hFE  as  a 
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function  of  IE  for  an  integration  interval  is  negligible. 


5.   MODEL  TESTS 

The  model  of  Fig.  12  has  been  exercised  in  a  single 
saturated  inverter  circuit,  in  an  emitter  follower  circuit, 
and  in  a  three  stage  d-c  coupled  amplifier  (5) .   The  saturated 
inverter  circuit  exhibits  the  behavior  of  the  model  under 
large  signal  application.   The  emitter  follower  circuit 
exhibits  the  behavior  of  the  model  under  small  signal  appli- 
cation.  The  three  stage  amplifier  exhibits  behavior  of  the 
model  under  both  small  and  large  signal  simultaneously.   The 
results  indicate  that  the  model  is  capable  of  predicting  in 
detail  both  the  steady-state  and  transient  response  of  the 
device . 

In  all  three  circuits  the  components  were  chosen  so  as 
not  to  obscure  the  device  performance.   Although  planar 
silicon  transistors  (Fairchild  1312)  were  used  in  the  experi- 
mental work,  the  model  has  more  general  application.   The 
systems  of  simultaneous  nonlinear  differential  equations 
resulting  from  these  circuits  have  been  programmed  on  the 
7090  in  conjunction  with  the  described  program  in  part  4  of 
this  paper. 
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5.1   CIRCUIT  CHARACTERIZATION 

The  circuits  used  for  these  studies  were  carefully  con- 
structed to  minimize  stray  inductance  and  capacitance.   The 
remaining  scray  capacitances  present  were  measured  in  the 
absence  of  the  transistor  and  added  to  respective  header  and 
land  capacitances  to  form  Cq-qq,    Cebs,  and  CCES.   The  remaining 
circuit  elements  were  measured  by  entirely  conventional  means. 
The  equivalent  resistance  of  the  pulse  generator  and  its 
termination  were  added  to  the  lumped  resistor  in  the  base 
lead  to  form  RQ.   A  carefully  calibrated  sampling  scope 
driving  an  X-Y  plotter  was  used  to  measure  the  .7  nanosecond 
ramp  input  pulse  and  the  output  voltage.   Shunt  capacitance 
of  the  scope  probes  were  taken  into  account  where  appropriate 
(5). 


5.2   SATURATED  INVERTER  CIRCUIT 

The  saturated  inverter  of  Fig.  A  was  driven  under  the 
conditions  of  the  resistance  RG  with  a  value  of  525  ohms,  vint 
was  a  pulse  with  the  peak  value  of  2.5  volts,  and  the  rise 
time  of  the  .7  nanosecond.   Fig.  18  programmed  on  the  7090 
computer  as  is  described.   The  system  of  nonlinear  differential 
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Fig.  18.   Nonlinear  transistor  model  in  saturating 
inverter  circuit. 


(Gg  4  GB^  "  GBVBP  +  (CCBS  +  cCEs)VB  "  CCBS^O 

=  GgVin(t)  (62) 

-GBVB  +  (GB  4.  GD  4  Gs)VBp  -  GSV0  4  (CE  •*  CD  +    Cc  4  Cs)VBp 

-(Cs  +  Cc  +  CCBS)V0  =  0  (63) 


-GsVBp  4-  (GL  +  GS)VQ  -  CCBSVB  -  (Cc  +  Cs)VBp 

+  (°CBS  4  CL  +  CC  +  CS^0  =  -hFEVBPGD  4  ECCGL 
•J-  hpEiGs (Vq  -  VBp) 


(64) 
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equations  describing  Fig.  18  are  given  in  equations  62  -  64. 

The  input  voltage  to  the  computer  was  described  as  shown 
in  Fig.  19  and  the  functional  relationships  representing  Fig. 
19  were  given  in  equations  65  -  71. 


Vin  =  VK  0<T  <TQ  (65) 

Vin  =  VK  *  VR(l-e"S(T-T0>)    TQ  s<  T  <  TR  (66) 

V±n  =  VL  TR^  T<TF  (67) 

Vin  =  VL  4  il-e^-^hVp      Tp<T<TQ  (68) 

Vin  =  VK  T  «TQ  ^69) 

Vr  -  VL 

X=   -m(yR  -  y\>)  (70) 

TR  "  T0 

1n/VL  ~  VF^ 

Tq  -  TF 


The  comparison  of  experimental  and  predicted  result  for 
saturated  inverter  Fig.  18  is  shown  in  Fig.  20.   The  experi- 
mental performance  of  the  device  agreed  with  the  predicted 
values  better  than  10  percent  for  turn-on  delay,  turn-on,  and 
turn-off.   Only  the  prediction  of  storage  time  did  not  agree 
with  experimental  results.   This  will  be  further  discussed  in 
section  6  of  this  paper. 
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Fig.  19.   This  form  of  the  input  is  simulated  on  the 
Computer. 
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T  —  4  nanosecond/division 


Fig.  20.   Comparison  of  experimental  and  predicted 
results  for  saturated  inverter. 


Solid  line  is  experimental  results. 

Dotted  line  is  predicted  results. 

V.   =2.5  volts  with  .7  nanosecond  rise  time. 

R    =525  ohms. 

y 
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5.3   EMITTER-FOLLOWER  CIRCUIT 

The  stability  of  the  emitter-follower  has  been  investi- 
gated by  many  authors.   It  has  been  pointed  out  that  the  input 
of  a  transistor  in  the  emitter-follower  configuration  appears, 
under  certain  conditions,  to  be  a  capacitance  shunted  with 
negative  resistance.   Therefore,  if  the  driving  source  to 
this  device  appears  to  be  inductive,  instability  and  oscilla- 
tion are  quite  likely  to  occur.   Fig.  21  presents  an  equiva- 
lent circuit  for  a  simple  emitter- follower .   The  device  model 
is  the  same  as  that. used  in  the  saturated  inverter  circuit 
with  the  exception  that  G   and  Cs  were  set  to  zero.   The 
deriving  source  contains  the  inductance  L  necessary  for 
instability.   This  inductance  L  in  practice  is  usually  caused 
by  lead  wire  inductance.   Fig.  22  shows  a  comparison  of 
theoretical  and  experimental  results  for  a  representative 
device  (transistor  X)  and  for  a  value  of  inductance.   The 
device  value  used  to  obtain  the  emitter-follower  results  were 
precisely  those  used  to  obtain  saturated  inverter  of  Fig.  20. 
This  result  demonstrates  that  this  model  need  not  be  res- 
tricted to  large  signal  applications  but  is  equally  useful 
in  representing  the  device  in  the  "linear"  region. 
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=  2.5  volts 


Fig.  21.   Nonlinear  transistor  model  in  emitter- 
follower  circuit. 
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Fig.  22.   Comparison  of  experimental  and  predicted 
results  for  emitter-follower. 


L   = 
Rn  = 


160  nh 
250  ohms 


r   =50  ohms. 


Vin(t) 


=  1  volt 


45 


5.4   THREE  STAGE  D-C  COUPLED  AMPLIFIER  CIRCUIT 

The  model  was  exercised  further  in  a  three  stage  d-c 
coupled  amplifier  of  Fig.  23.   The  device  model  was  the  same 
as  that  used  in  the  saturated  inverter  circuit.   The  circuit 
used  for  this  study  was  carefully  constructed  to  minimize 
stray  inductance.   Stray  inductances  due  to  the  circuit-wire 
were  found  to  be  negligible.   No  inductance  was  programmed 
as  it  was  done  for  emitter-follower.   Fig.  24  shows  a  compari- 
son of  the  theoretical  and  experimental  results  for  a  repre- 
sentative device  (transistor  X) .   The  result  indicates  that 
the  use  of  a  nonlinear  model  makes  it  possible  to  obtain  both 
the  steady-state  and  transient  solution  simultaneously. 

In  the  analysis,  T. ,  T2  and  T3  are  represented  by  the 
nonlinear  transistor  model  of  Fig.  12. 
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Fig.  23.   Three-stage  D-C  Coupled  Amplifier 
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Fig.  24.   Comparison  of  predicted  and  experimental 
results  for  the  circuit  of  Fig."  23. 


V   =  5  mV  pulse,  500   sec.  in  duration, 
in 
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6.   MODIFICATION  OF  THE  MODEL 

The  extrinsic  model  of  Fig.  12  is  capable  of  predicting 
the  turn-on,  turn-on  delay,  turn-off  delay,  and  d-c  characteris- 
tics of  the  transistor  in  detail.   However,  it  fails  to  pre- 
dict the  turn-off  (storage  time)  in  detail.   This  suggests 
that  the  extrinsic  model  should  he  further  completed.   Fig.  25 
is  a  plot  of  collector  deplation  capacitance  versus  junction 
voltage.   The  equation  (49)  is  used  to  describe  this  phenomena. 
The  result  indicates  that  the  equation  (49)  is  only  capable 
of  the  predicting  of  the  depletion  layer  capacitance  for 
reversed  biased  condition.   This  indicates  that  the  functional 
relationship  used  to  describe  the  deplation  junction  capaci- 
tance was  inadequate  and  failed  to  predict  forward  biased 
characteristic.   This  failure  used  to  be  blamed  on  the  1) 
difficulty  of  reproducing  the  forward  biased  data  characteristic 
of  deplation  capacitance,  2)  the  diffusion  capacitance.   A 
careful  capacitance  measurement  with  a-c  signal  as  small  as 
20  microvolts  has  made  it  possible  to  reproduce  the  forward 
biased  data.   The  diffusion  capacitance  values  were  calculated 
and  they  indicated  that  for  emitter  junction  the  calculated 
value  of  the  diffusion  capacitance  were  very  close  to  the 


48 


observed  values.   However,  for  collector  junction  the  diffu- 
sion capacitance  was  by  order  of  magnitude  smaller  than 
observed  value.   As  a  result,  a  new  functional  relationship 
for  describing  deplation  layer  capacitance  is  being  adopted. 
The  new  functional  relationship  is  capable  of  predicting 
depletion  layer  capacitance  for  the  forward  biased  condition 
as  well  as  reversed  biased  condition  within  five  percent  of 
the  measured  value.   The  results  are  plotted  in  Fig.  25  and 
the  functional  relationship  is  given  by  equation  (7  2) . 


C   =     CC1        +      CFC1  ,72) 

cc   iv^c  -  vCB)iU   (v^c  -  vCB)^       (72) 

Where  Cq^,    m,  CpC-|_  and  p  are  found  empirically. 


Although  this  new  functional  relationship  is  not  in  exact 
exponential  form,  it  closely  resembles  the  work  done  by  C.  T. 
Sah  and  will  predict  the  same  result  at  the  limits.   C.  T. 
Sah  attributes  this  effect  to  free  carriers  in  the  transistor 
region  (8) . 

The  new  functional  relationship  makes  it  possible  to 
predict  storage  time  within  a  few  percent.   The  results  are 
plotted  in  Fig.  26.   The  extrinsic  model  is  called  complete 
model  when  substituting  equation  (72)  for  equation  (49)  . 
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C    =      CC1       4      CFC1 
CB    (V^c  "  VC3)ra    (V^C  ~  VCB)r 
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©©©  CEB  Experimental 
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(V^  -  V)  Volts 


t. 4.  5'.  bl  j 


Fig.  25.   Collector  depletion  layer  capacitance  v.s 
voltage  across  the  junction. 


50 


T  —  4  nanosecond/Division 


Fig.  26.   Comparison  of  predicted  and  experimental 
results  for  the  circuit,  Fig.  18. 


Solid  line   =  experimental  results 
Dotted  line  =  predicted  results. 
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7.   SENSITIVITY  OF  PREDICTED  RESULT  TO  DEVICE 
AND  CIRCUIT  PARAMETERS 


Every  measurement  made  is  subject  to  a  number  of  errors. 
In  order  to  determine  the  influence  of  measurement  errors  on 
the  predicted  results,  nominal  values  for  each  parameter  were 
used  in  the  inverter  program  (5) .   These  results  were  compared 
to  the  results  obtained  when  each  nominal  value  was  incremented 
independently  by  5  percent.   The  sensitivity  results  indicate 
that  Cp  and  the  parameters  that  determine  base  current  V.  , 
RG  ■*■  RB'  an<^   were  by  far  the  most  important  parameters  for 
this  circuit.   The  circuit  response  is  relatively  insensitive 
to  parameters  such  as  hFE  and  T -^.      The  results  of  the  sensiti- 
vity study  make  it  possible  to  determine  which  device  parameters 
need  only  be  specified  nominally,  thus  reducing  drastically 
the  number  of  parameters  to  be  measured. 

The  program  from  which  these  results  were  obtained  has 
facilities  for  emitter  current  dependence  of  hFE  and  R„. 
However,  all  the  results  presented  assume  constant  values  for 
hFE  and  Rg.   The  predicted  result  from  nonlinear  hFE  and  Rs 
agree  to  the  third  place  with  the  linear  results.   If  a  circuit 
with  a  saturation  current  approaching  ma  for  the  transistor  X 
had  been  chosen,  it  would  be  anticipated  that  results  would 
be  a  great  deal  different. 
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8.   CONCLUSION 

The  results  shown  in  this  paper  demonstrate  the  feasi- 
bility of  using  a  nonlinear,  equivalent  circuit  model  to 
represent  a  transistor  in  large  signal,  small  signal  and  both 
large  signal  and  small  signal  operation  simultaneously. 
Although  the  range  of  devices  and  circuit  conditions  is  quite 
limited  in  this  study,  it  is  significant  to  note  that  none 
of  the  devices  strayed  more  than  ten  percent  from  their  pre- 
dicted behavior,  and  the  majority  of  this  is  attributed  to 
the  great  dependence  the  circuit  has  to  V-  n(j-\  . 

The  model  does  not  require  any  advanced  programming  and 
machine  running  time  is  not  lengthy.   Machine  running  time 
for  the  results  of  saturated  inverter  circuit  was  forty 
seconds,  including  read  and  write  tape. 

The  model  does  an  excellent  job  in  predicting  transient 
time.   The  model  is  not  a  cumbersome  model.   It  does  not 
introduce  any  new  circuit  elements.   It  is  compatible  with 
the  other  circuit  elements  to  which  it  is  to  be  connected 
and  it  does  not  require  a  circuit  designer  to  be  accustomed 
to  new  circuit  elements.   The  relationship  to  device  physics 
is  relatively  clear,  and  the  model  is  capable  of  providing 
device  people  with  information  about  their  products. 
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The  linear  form  of  the  model  has  found  wide  acceptance 
in  the  design  of  small  signal  amplifiers.   The  introduction 
of  nonlinear  elements  intensifies  the  value  of  this  model  for 
small  signal  design  since  it  permits  simultaneous  a-c  and 
d-c  analysis. 

Parameters  of  the  model  can  be  evaluated  easily  with 
conventional  techniques.   Circuit  elements  and  device  para- 
meters may  be  optomized  long  before  the  devices  are  available. 
This  approach  is  particularly  useful  in  the  design  of  inte- 
grated circuit  where  the  trimming  of  a  circuit  is  impractical. 
The  model  can  be  used  as  the  basis  of  an  integrated  specifi- 
cation which  would  include  large  signal  and  small  signal 
application. 

The  model  and  the  paragram  can  be  combined  to  a  subroutine 
paragram.   The  subroutine  can  be  used  instead  of  transistor 
in  the  circuit,  and  the  performance  of  the  circuit  can  be 
predicted  easily. 

In  short,  this  model  is  a  practical  model. 
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APPENDIX  A 

Transistor  X  is  a  NPN  planar  silicon  transistor  manu- 
factured by  Fair child  and  was  designated  type  2B8. 

The  breakdown  voltage  of  either  junction  is  about  9 
volts.   The  hFE  of  the  transistor  is  about  75,  the  detail 
hFE  characteristic  is  given  in  Fig.  16  of  this  paper.   Other 
characteristics  of  this  transistor  are  given  in  the  device 
characteristics  of  this  paper. 

Other  types  of  transistors  were  used  in  the  experimental 
work,  which  all  have  exhibited  similar  results  as  transistor 
X.   However,  all  the  data  and  results  represented  in  this 
paper  are  of  transistor  X. 
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APPENDIX  B 

The  listing  of  the  program  used  to  predict  the  theore- 
tical value  is  given  in  this  Appendix.  The  data  and  pro- 
gramming of  this  Appendix  is  a  portion  of  a  published  paper 
by  the  International  Business  Machines  Corporation  which  is 
a  supplement  to  reference  (X).  The  paper  identification  is 
as  follows. 

Identification 

Nonlinear  Transistor  model  for  the  prediction  of 
circuit  transient  response 

R.  J.  V7ilfinger/P.  Kiankhooy-Fard/S.  C.  Plumb 

Program  Description  -  A.  Selby,  Scientific  Computa- 
tion, International  Business  Machines  Corporation, 
Components  Division,  Fishkill,  New  York 
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*«*#>i-*##***-*tt#*##i*#-  LISTING  i********************** 

C      TRANSIST  V  CF  ADDITION 

C 

C      VERSION  WITH  HYPERBOLIC  FORM  OF  CF 

C 

DIMENSION  A(80)9  A I R ( 50 ) »  A I H ( 50 )  9  RB(50>»  HFE(50)»  V(4)»  VP(4)» 
S  VC(4) »  TBL( 16) 

EQUIVALENCE  (V(l)»  V8P)v  (V(2)9  V0)»  (V(3)»  VB)»  ( VP ( 1 ) »  DVBP ) » 
S   (VP(2)»  DVO)»  (VP(3)»  DVB)  »  (V(4)»  VG )  »  (DVG»  VP(4)) 
READ  INPUT  TAPE  2»  5»  A 

5  FORMAT  (20A4) 

1  READ  INPUT  TAPE  2»  2s  IVSN»  ERR 

2  FORMAT  (I5»  5Xs>  E5*0) 

READ  INPUT  TAPE  2»  3»  VG » VB 9 VBP » VO9DT »TMAX 

READ  INPUT  TAPE  2»  3»  VK»  VL&  VR»  TO?  TR 

READ  INPUT  TAPE  2»  3»  VF9  TF9  TQ 

READ  INPUT  TAPE  2s  3*  RG»  CBES»  CCBS?  RL »  CL»  ECC 

READ  INPUT  TAPE  2»  3s>   CCES  »R2  ?C2  9R3  »R  »C 

READ  INPUT  TAPE  2»  3?  AISE»  GAM9  Q»  P»  TB»  VPHI 

READ  INPUT  TAPE  2*    3»  CE1»  EM»  CC1»  EN,AA  »  BB 

READ  INPUT  TAPE  2*  3»AISC  »  AMU  »  VPHC».  TS»  HFE I  .HFX 
3.    FORMAT  (6E12o4) 

READ  INPUT  TAPE  2»  4»  Ni  (AIH(I)»  HFE(I)»  I  =  1»N) 

READ  INPUT  TAPE  2*  4?  M»  (AIR(I)t  RB (  I  )  9  I  ■  1>M) 
4     FORMAT  (15/  (2E12.4)) 

CALL  GIRL 

CALL  TITLE 

WRITE  OUTPUT  TAPE  3»  6»  A(l)»  VB»  A(2)»  VBP»  A(3)»  V0»  A(4)t  DT» 
$A( 5) 9  TMAX 

WRITE  OUTPUT  TAPE  39  69  A(6)»  VK9  A(7)»  VL9  A(8)»  VR»  A(9)»  TO9  A( 
S  1 0  )  9  T  R 

WRITE  OUTPUT  TAPE  39  69  A(ll>9  VF9  A{12)»  TF9  A(13)»  TQ9 A ( 55 ) 9CCES 

WRITE  OUTPUT  TAPE  3»  69  A(14)»  RG9  A(15)»  CBES9  A(16)t  CCBS9  A(17) 
$9  RL9  A( 18) 9  CL9  A( 19)  9  ECC 

WRITE  OUTPUT  TAPE  39  69  A(20)»  AISE9  A(21)»  GAMt  A(22)»  Q»  A(23)9 
$  P9  A(24)9  AISC9  A(25)9  AMU 

WRITE  OUTPUT  TAPE  39  69  A(26)»  VPHI»  A(27).  CElt  A(28)»  EM.  A(29)t 
S  CC19  A(30)9  EN9  A(31)9  HFEI 

WRITE  OUTPUT  TAPE  3?  6»  A(53)9  VG9  A(56)»  R.  A(57)9  Ct  A(58)»  R2f 
$  A( 59) 9  C29  A (60) 9  R3 

GG  =  l./RG 

GL  =  l./RL 

G  =  lo/R 

G2  =  1./R2 

G3  =  1./R3 

ALPHA  =  LOGF((VL  -  VF ) / ( VK  -  VF))/(TO  -  TF) 

DELTA  =  -LOGF((VR  -  VL)/(VR  -  VK))/(TR  -  TO) 

WRITE  OUTPUT  TAPE  39  6»'A(32}»  TS»  A(33>9  TB9  A(51)t  ALPHAt  A(52)» 
$  DELTA?  A(61)»  AA9  A(62)»  BB 

WRITE  OUTPUT  TAPE  39  69  A(64)9VPHC 

6  FORMAT  (IX9  6(A49  E12.5.  4X ) ) 
T  =  0, 

rc  =0. 
I  =  0 
NN  =  4 


7.  ; 

■ .    .    ■ 

—  -  -  - 

2  j 
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•V 


ro     <  - 


. 

.   - 

.  ,"     -     , 

. 

! '  s    VF    + 

1C 

.   1 

M' )     i  1 0  s     ■           10 

F     .--..'..'.,,      '  ,     -    TF)  ] 

••C      -•       /G    ■      :CC)    oL    +    DV0*CL    +GG*(V8    -    VIN) 

-      •     .       ■  .  •        20*130 

L  2C       -,:-••-. 
C    16 
D       >       14G    -    =    2»i 

:  C     i  .  I 
CONTINUE 

"•'■'•  ■  . 

..'    ro  i6C 

-  ■-.'    -  :     .   ■■  !  AIHi  1-1 )  )  /UIH(  I  )    -   AIH( 1-1 ) )*(HFE( I )    - 

i      HFEII-Iil 
30     "( 
6         iF(AI£-AIR<1.)        I  / ■      3      •  .•    6 0 
tR(  1 

:    i 
l^o     r>«      90  i  =  i 

IF    ( AIE    -  I) i  ■     >0G»19( 

I  N  U  E 

i     ;    .    f 
■ 

-1)     h     (AI  IRi  1-1)  )/    MR(  I  )     -   AIR( 1-1) )*(R8( I)    -    RBU 


:       '-■:-' 

■        .  .      - ,      p   ■    vc 

21  .  ■     -        [SE       -  •-.  <  Hi    :+l,  }*V8P) 

ISE 

=    C  E  I  •       ■  -     .  E 

c  "   -       / ,  •       ;    ■     '  0         ,-,b 

■--    CC1/       0    -        1       +    VP  K 

[Y-l. ) )/(     HFEI  +    1 .. )*(VBP- -VO) ) 

c  -   .   s     -.  iu     &  isc     n/<   hfei  +] ,       i  :f 

I  '.    =       »     *B 
;A     ■    (  VB    -   V3  ?  5  -;.-GE     *■    V8    ^2    +    [  VG    -    V I N  5  *     3 

•     /P,P*(G8     r      ,n    ■     >>c.  .     -    V  ! » "  -    -    < 
X.     i    VO*(GL    -:-    GS3  V8P*GS    +  ••;•'.-."'.     ,;;    -    £CC*GL    -    r,,rti^GS*(VbP    - 

)J  VO*G"$ 

X&"CBES+CCB$*C2 

-  ce  -    cr   -r  cc   • 

=   cc  -.-   cs 

X^     k    CCBc     +■    CL    +    CC     i"    C       -     ' 

>V0    -    -i     "■  •  .  .   '     •      '        II/XG)    •;•    J      '■>;'-  2     »■    XD*(-XI**2/XG    r    X J )  ) 
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An  approach  to  transistor  transient  and  steady  state  pre- 
diction for  large  signal,  small  signal,  and  both  large  signal 
and  small  signal  simultaneously  is  described  which  uses  a  non- 
linear equivalent  circuit  model.   The  model  is  based  on  the 
Eber  and  Moll's  equations  (intrinsic  model).   The  model. is 
modified  by  introducing  extrinsic  elements.   The  model  is 
further  completed  by  further  investigation  of  collector  junc- 
tion under  forward  biased  condition.   A  comparison  of  experi- 
mental and  predicted  results  are  presented  for  a  saturated 
inverter  circuit,  an  emitter-follower,  and  a  three  stage  d-c 
coupled  amplifier.   The  method  of  device  characterization  is 
outlined  in  detail.   A  sensitivity  study  is  performed,  and  the 
relative  sensitivity  of  the  device  parameter  is  pointed  out. 
Theoretical  results  are  obtained  from  numerics  solution  of  a 
system  of  nonlinear  differential  equations  describing  the 
model  and  circuit  diagram.   The  nature  of  the  errors  in  the 
program  are  pointed  out.   The  capability  of  the  model  to  linear 
amplifier  design  is  discussed  and  an  integrated  switching  and 
linear  amplifier  specification  is  proposed.   A  subroutine  pro- 
gram is  proposed,  which  will  enable  the  circuit  designer  to 
handle  the  transistor  as  one  would  have  handled  such  circuit 
elements  as  the  resistor,  or  the  capacitor  in  a  circuit 
without  much  knowledge  about  the  material  of  the  resistor  or 
the  capacitor. 


